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EXACT SOLUTIONS OF NONSTATIONARY 
SHRODINGER EQUATIONS and GEOMETRIC PHASE 
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^ I Abstract. The procedure of solving nonstationary Schrodinger equations in the 

exact analytic form is elaborated on the basis of exactly solvable stationary models. 
The exact solutions are employed to study the nonadiabatic geometric phase at 
cychc evolution of quantum systems. 
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1 Introduction 



> 

f^ ' Consider the time-dependent Schrodinger equation 

§ : m?m2L . Hitmt) > (1) 

0> I the exact solution of which is not known in general. The purpose of this paper is 

^ ' to give the procedure of obtaining a wide class of Hamiltonians for which the exact 

'P' solutions of ([I|) can be found. 

'S . Introduce a unitary transformation 5^(t) = <S^^{t) 

a ■ 

gn: \^{t)>=S{tMt)> (2) 



such that in the new equation 

~ dt 



S : a^^#^ = sm > (3) 



the Hamiltonian 

H = S^{t)H{t)S{t) - ihS\t)S{t) (4) 

does not depend on time, i.e., {d/dt)H = ||l|, |^. From this, we obtain the equation 
for the unitary transformation S(t) 

SW{t)S + S\dH{t)/dt)S + S^H{t)S - ihS^S - ihS^S = 0, (5) 

where S = {d/dt)S. 

From (^) for t = it follows 

H = S^{0)HoS{0) - ihS\0)S{0). (6) 
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Taking account of this equation in (^) we obtain 

H{t)=S{t)S\0)HoS{0)S\t), (7) 

where Ho = H{t = 0). 

Let us define the solution to equation (|^). On the one hand, it is 

r){t) = ex.-p{-iHt/h)r]{0), (8) 

on the ather hand, it is 

ri{t) = exp{-iHt/n)<!>{S). 

Therefore, 77(0) = $(^) is the function of the initial state of (^ and is one of the 
eigenfunctions of H 

H<^{S)=S<t>{S). (9) 

Now the solution ([^) is reduced to the solutions of equations (^ and (|^). The 
solution of the equation (|) becomes the eigenequation problem (||) for the time- 
independent Hamiltonian H. But solving (^) is in general a very complex problem. 
On the other hand, it is possible to assume the concrete forms of S{t) and seek the 
Hamiltonian H{t) that satisfies (^. Then, to obtain exact solutions of (|I]) we, in 
fact, need an exactly solvable time-independent equation (^. After this, by using 
the concrete forms of S{t), we can construct a wide class of exactly solvable time- 
dependent H{t) 0]. 

In terms of the time-evolution operator U{t) = U(t,0) defined by, \E'(t) = 
W(t)\l/(0), the Schrodinger equation (|1|) is written as 

inU{t) = H{t)U{t). (10) 

It is not difficult to derive a very important relationship between the operators U{t) 
and S{t) 

U{t) = S{t) exp{-iHt/h)S^{0) (11) 

because if H and S{t) are known, U{t) is known, too. 

2 Construction of the Hamiltonian of the non- 
stationary Schrodinger equations 

Consider an example of reconstruction of a time- dependent Hamiltonian 

H{t) = {pl + q{x))i + Bit,x)-j (12) 

for which equation (|l|) is exactly solved. Here px = —ihVx is the operator of 
momentum, q{x) is a stationary potential and B(t, x) is a time- varying potential that 
can be considered as an analog of the inhomogeneous magnetic field, j = {jx,jy,jz) 
is the spin operator and ji = (2j + 1) x (2j + 1) are matrices. 



Let the model Hamiltonian H{x) of the stationary problem correspond to two- 
channel one 

fiVnix)-V22ix))/2 Vuix) \ 

^ \ V^2i(x) -iV^^ix)-V22ix))/2J■ ^'^^ 

It can be rewritten as 

H{x) = {pl + q{x))i + B{x)-CT (14) 

sin 6{x) 



{p^ + q{x))I + n{x) , ■ ni \ at \ 

^ Vsm6'(x) — cos6'(x) 

with the evident notation: g(x) = {Vu^x) + V22(a;))/2, / is the unit matrix, cr = 
((Ti, (72, 0-3), B(a;) is expressed in terms of Vij{x) and can be defined as 

B(x) = n(x)(sin^(x),0,cos^(x)), (15) 



n(x) = l/2^(Vn(x) - V22{x)f + AVi2{x)\ 

""^^") = T](^' ^°^'(") = Wx) ■ 

By means of a canonical unitary transformation S{t\ one can obtain a time-dependent 
Hamiltonian Hii) in an explicit form from the stationary Hamiltonian (0). If Sif) 
is chosen as an operator of rotation around 2;-axis 

S{f) = exp(— iascut), (16) 

then in accordance with (^) the following relationship between the stationary H and 
nonstationary H[t) Hamiltonians is obtained 

H{t) = S{t)H{x)S^{t) + uja3. (17) 

This time-dependent Hamiltonian with an allowance for (|1^ can be rewritten as 

cos 6{x) + uj/Q{x) sin 6{x) exp{—2iujt) 



H{t) = {{p^ + q{x))I + ^{x) y ■ at \ /,o- j.\ / /]/ \ , ifM \\ } ■ 
\j \\i-x IV // V ^ \sm9{x) exp{+2iujt) —{cos9{x)+uj/il{x))J 

This may correspond to the problem of two-level atom for the spin-1/2 representation 
of j. Obviously, it corresponds to the Hamiltonian (|T2|) with the nonstationary field 

B(t,a;) = r2(a;)(sin6'(a;) cos(2co't), sin^(x) sin(2co't), cos^(x)), (19) 

the components of which are expressed in terms of the components of the time- 
independent field 

Q{x) = Q{x){l + ^ + 2^cos^»)'/', (20) 

^. , r2(a;) sin^(s) r2(x) cos^(2;) u 

sm6'(x) = — r^; , cost^(a;) 



n{x) ' ' ' n{x) n{x)' 



As a result of the unitary transformation S{t), the -82(^5 x)- component of the vector 
B(t, x) arises. The expression for the initial Hamiltonian Ho{x) is directly obtained 
from fllTD or from fO) at t = 



Ho{x) = H{x) + uoa^. (21) 

This relation is conveniently rewritten in the form 

Ho{x) = iPl + q{x))i + Bo{x) ■ cr, (22) 



where the potential Bo(a;) coincides with B(t, x) defined in (|T9D at the initial moment 
t = 

Bo{x) = fi(x)(sin^(x),O,cos0(x)). (23) 

In accordance with (|1^) and (|21|) and with allowance of ([16|) the time- varying Hamil- 
tonian H{t) is written in terms of the initial Hamiltonian as 

H{t, x) = exp{-ia3Ujt)Ho{x) expl+ia^ut) = {pl + q{x))i + B(t, x) ■ cr. (24) 

It means that the time-dependent part of the reconstructed potential is in precession 
along the 2;-axis with the frequency 2u. Thus, one has been built the time-dependent 
Hamiltonian 

Hit) ^ (p^ + ^^iM±M^)/ 

/ viiM_v^aM + ^ V,,{x) eM-2^cot) \ ..,^ 

+ [v2^{x)eM+2^u^t) -(^^iM_lM£) + ^) j ^'^) 

that admits exact solutions of the nonstationary Schrodinger equations (Q) if the 
corresponding canonical stationary equation (H) with the Hamiltonian (0) is exactly 
solved 

\'^{t,x) >= exp{-ia3Lot) exp{-iHt/h)\(^{x) > . (26) 

Indeed, the eigenvalues S of the stationary equation with the Hamiltonian H are 
known because they are initial data for the inverse scattering problem. By using 
the algebraic procedure of the Bargmann-Darboux transformations, the potentials 
Vaa'{x) in ([l3|) and solutions |$(x) > are reconstructed in the exact analytic form 
from the known spectral data. 

Let V{x) be a Bargmann potential matrix for which the system of equations (|^) 
has exact solutions. For example, matrix elements of the transparent potential and 
pertinent Jost solutions are Q 

K„'(a:) = 2^Y1 ^M-<xh:Pux{xha' exp(-4x), (27) 

Paa'i^^x) = exp{±ikaX)5aa' (28) 

- ^7>xp(-<x)p-/(x) / ^^,exp{-{Kl,±tk^,)x')dx\ 

uX -^^ 



where 

The normahzed solution for the vector function of the bound state is immediately 
expressed through the matrix Jost solutions taken at the energy of the bound state 

m 

\^a,{£v,x) >= J2P»a'{£u,x)'j'^,, a, a' = 1,2. 

a' 

Here, the indices u and A correspond to the bound states characterized by the 
energies Si, and normalizing matrices | 7q >< 7^/ 1 and the channel indices are labeled 
as a, a'. It is worth noting that the class of exactly solvable stationary problems 
is very wide. Algebraic Bargmann-Darboux transformations give multiple examples 
of such problems that can be used to obtain the exact solution of the Schrodinger 
equations for a large class of explicit time-dependent Hamiltonians. 

3 The geometric phase 

By means of the unitary transformations S{x) and S{x) one can annihilate the 
Bi{x)- and i?i(x)-components of the fields B(x) and Bo(x), respectively. It is 
realized by means of the S'f/(2)-operators 

S{x) = exp{—i9{x)J2) and S{x) = exp{—i9{x)J2). (29) 

Since the following relations are valid: 

B(a;) ■ <T = fl{x)exp{—i9{x)J2)d'3exp{i9{x)J2), (30) 

Bo(x) ■ a- = Q{x) exp{-i9{x)J2)cr3exp{i9{x)J2), (31) 

under these transformations S{x) and S{x), the systems of equations with the Hamil- 
tonians ( [T^ ) and ( P^ are reduced to the gauge-type equations [Q 

[-{Vj + A{x)Y + q{x)i + n{x) ■ as - P^]^\x) = 0, (32) 

[-{Vj + A{x)y + qix)I + Q{x) ■ a, - P']%{x) = 0, (33) 

P = diagipa) 

for the new functions $'(a;) and '^'o{x) 

$'(x) =exp(-z^(x)j2)"I'(x), ^^(3;) = exp(-z6'(a;)J2)^o(a:). (34) 

The matrix elements of the effective vector potentials A{x) and A{x) in ( |32D and 
(|33|) are generated by the procedure of the gauge transformation (pOJ) 

A{x) = exp{i9{x)J2){d/dx) exp{—i9{x)J2), (35) 

A{x) = exp{i9{x)J2){d/dx) exp{—i9{x)J2). 



Note, the elements of S{x) 

^ ^ ^v V yj'^/ \-sm6{x)/2 cos6{x)/2j 

are defined by the elements of the known potential matrix Vaa'{x) (|l5|). Taking into 
account (|^) for the transformation operator S(t) in ([TT|), we get 

U{t) = exp{—id'3Ujt)exp{—iHt/h). 

Then, the solution (^) of (^ with the Hamiltonian (^) can be rewritten as 

The eigenstates and eigenvalues of H{x) are \^a{£u) > and S^ = £^I, a = 1,2. 
After one period T = n/uj, the solution is 

\^a{T,S,) >= expiTin) exp{-tS,T/n)\^^{0,S,) > . (36) 

and the total phase is 

5, = i±7r + £^T/h). (37) 

Next, to find the dynamical phase we calculate the expectation value of H{t) by 
using the relations (0) and (||) or (|2! 



=< <!?^i£,)\Hix)\^a{£u) > + < ^ai£u)\uja3\<!?^{£,) >= £, + ud^. (38) 

By virtue of ([3^ ) for the eigenf unctions of H{x) the spin-expectation value a^ is 

written as 

dl =< ^a{t,£u)\cr3\'^a{t,£u) > = < $a(^.) |a3|$a(^.) > (39) 

=< (^'^{£^)\exp{ie{x)j2)a3exp{-ie{x)j2\^'^{£u) > 

=< <^',{£,)\ - sm{e{x))ai + cos{e{x))as\K{^,) > . (40) 

Then, for the dynamical phase 6"^ we get the expression 

T 
St = J <^a{t)\H{t)\^^{t)>,dt 


= £,T/h + TT < <^'^{£,)\ - sin(^(x))ai + cos{e{x))a3W^{£,) > . (41) 



Removing the dynamical phase (|41| ) from the total phase ([37| ) we get the relation 
for the nonadiabatic geometric phase 51 

5l = {5,-5t) = 7r<$:,(£,)|sin(^(x))ai + [±l-cos(^(x))a3)]|$L(^-)> 

= vr(±l-^3). (42) 



One can easily see that the geometric phase is determined by the spin-expectation 
value o"3 along the rotating axis. For the first time this result was obtained by 
Wang in |l[] for a spin particle in a rotating magnetic field. As distinct from |l|] this 
approach is based on the exact solutions |$a(x) > of the body-fixed Hamiltonian 
(JT^). In the case when the dependence on the space variable is absent, our problem is 
simplified and goes to the problem for the spin-1/2 particle in the rotating magnetic 
field with the cranking Hamiltonian 

H(t) = exp{—icr-^ut)Hoexp{+ia3UJt) = B(t) ■ cr, (43) 

B(t) = r2(sin^cos(2ti;t), sin^sin(2co't), cos6) 



instead of the inhomogeneous H{t,x) in (p^) or ([T8|). It is evident from (^) that 
the Routhian (the body- fixed Hamiltonian H{x), ([I^ ) ) turns into 

H = Ho-uja3 = B-cT (44) 

with the renormalized magnetic field 

B = 0(sin^,0,cos^). 

The initial Hamiltonian Ho of this model is 

Ho = Bo-cr, Bo = fi(sin^,0,cos^). (45) 

The total phase is determined as in (^) and the dynamical phase can be obtained 
from ( p8D taking account of ( ^3]) for H{t) and with the spin alignment a^ determined 

as 

^^ =< ^mit)\a3\^mit) > = < ^jm\cr3\^jm > (46) 

=< jm\ exp{ieJ2)d-3exp{-i9J2\jm > 

=< jm\ — sin(^)ai + cos{9)d-3\jm >= mcos{9), (47) 

where \jm > are eigenfunctions of a^, the channel index is m instead of a. Thus, 
the dynamical phase is 

^m = £mT/h + mn cos{6) 
and the geometric phase is 

(5^ = m7r(l - cos(^)) (48) 

In the adiabatic limit, when u/Q —>■ 0, then Q —>■ Cl, cos6 —>■ cos6 and the geometric 
phase (^Sp becomes Berry's phase ^ 

5^ = m7r(l — cos^). 

Note, the relation (^) corresponds to the nonadiabatic Aharonov-Anandan geo- 
metric phase 

T T 

SB = fdt< '^'{t)\{d/dt)-^'{t) >= fdt< <l>„(x)|5"f(t)5(t)|$„(x) >= Tcu^s = 71^3.(49) 



We have used the stationary two-channel exactly solvable models and a quite simple 
time-dependent unitary operator S(t) of rotation to construct the time- dependent 
Hamiltonians and the corresponding solutions of the nonstationary Schrodinger 
equations in a closed analytic form. It is evident, how the suggested approach is gen- 
eralized for another choice of stationary models and time-dependent iS(t)-operators. 



4 Conclusion 

The procedure is given of constructing time-dependent Hamiltonians for which non- 
stationary Schrodinger equations have exact solutions. Time-dependent unitary 
transformations and exactly solvable models for stationary Schrodinger equations 
are used to obtain a large class of explicit time- dependent Hamiltonians. The two- 
channel problem has been considered in detail, and exact solutions are used to 
calculate the nonadiabatic geometric phase. 
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